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In these lectures I discuss the pos-ible consequences of forming cosmic 
strings and superconducting cosmi< strings in the early universe. Lecture 
1 describes the group theoretic reasons for and the field theoretic reasons 
why cosmic strings can form in pontaneously broken gauge theories. 
Lecture 2 discusses the accretion of matter onto string loops, emphasising 
the scenario with a cold dark matter dominated universe. In lecture 3 
superconducting cosmic strings are discussed, as is a mechanism which 
leads to the formation of structure from such strings. 
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1 Lecture 1: Cosmic Strings and Phase Transi- 
tions in the Early Universe 

Recently, astronomers have been syster latically probing the very large scale struc- 
ture of our universe, and have made observations that theorists are now trying to 
explain. For example, what causal mec lanism could produce the primordial energy 
density perturbations which are thought of as necessary to seed galaxies and clusters 
of galaxies? How could the inhomogene ity represented by this large scale structure 
and the galaxies be reconciled with the observed smoothness of the microwave back- 
ground radiation? The interface of paricle physics and cosmology has provided us 
with one of the most intriguing possible s flutions to this problem. The idea that phys- 
ical processes occurring just secoi ds after the initial big bang should directly 

determine the structures being observec some 15 billion years later is staggering in 
the extreme. Yet this is just what two re 'ent particle physics ideas, used to determine 
the initial density fluctuations, suggest. The first is based on quantum fluctuations 
which arise in Inflationary universe mot els 1,2 . I will not discuss these models, but 
strongly advise anyone interested in cosmology to refer to the excellent review arti- 
cles already written on the subject. A1 credits' article 2 gives a very clear physical 
picture for the reasons behind inflation The second is based on cosmic strings 3 . 
and a review of their properties is the g< ai of these lectures. 

In this lecture I will discuss how the strings form and how they evolve in an ex- 
panding universe. In lecture 2, I will describe how initial density perturbations grow 
around long strings and string loops emphasising the Cold Dark Matter (CDM) sce- 
nario. In my final lecture I will discuss i very exciting model due to Witten[4;, who 
realised strings may be superconducting in that they can carry persistent currents. 
This model was later used to account foi the formation of the large scale structure in 
our universe, although rather than acting as seeds around which matter can accrete, 
these superconducting cosmic strings a :t as seeds for explosions as they emit elec- 
tromagnetic radiation, forming bubbles of plasma, on the shells of which the galaxies 
form and fragment [5]. We will see how -he scenario is restricted by the dynamics of 
such loops, indeed to the extent that it a ppears the scenario works only for a limited 
range of coupling constants. 
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In motivating the role of particle physics in modern cosmology, I refer the reader 
to some of the excellent review articles written on the subject[6]. The hot big bang 
theory of the early universe successfully predicts the Hubble expansion, the microwave 
background radiation and the light element abundances. It appears to fit in nicely 
with ideas of particle physics where, as the energy is increased so is the degree of 
symmetry used to describe the particle interactions. At high enough energies (i.e 
the very early universe when the temperature was very hot) we find the universe in 
the state of maximum symmetry. This is spontaneously broken as the universe ex- 
pands and cools through some critical temperature. However what is required is some 
source of the perturbations essential to produce the structures we see today. Amongst 
the more unusual large scale features are giant ‘filaments 1 (linear overdense regions 
in the galaxy distribution, about 200ft fg 1 Megaparsecs long and lOftjgMpc across 7 , 
large ‘voids 1 , empty of bright galaxies. 120fir r) 1 Mp C in diameter, and galaxies lying on 
the surface of 'bubbles' dO-COftjgMpc across S'. (ft 5G is Hubbles constant in units of 
50 kms~ l ). An Abell cluster is defined to be a region smaller than 3ftjgMpc in radius 
containing more than 50 bright galaxies. These clusters appear to be clustered on 
scales of lOOh^/Mpc with a mean separation of llOft^Mpc. 

Now it appears that gravity alone could not have moved galaxies and led to such 
large scale structure since the big bang. Turok’9: explains this in a succinct argument 
which I will follow here. Peculiar velocities (velocities relative to the Hubble flow) 
grow as t j in an expanding universe. As we shall see in the next lecture, in the linear 
regime, ^ < 1, there is a precise relation, Sr = Hq 1 Sv w'here Sr is the peculiar dis- 
placement, Sv is the peculiar velocity and the Hubble radius Hq 1 , (characterising the 
expansion rate of the universe) = fiOOOh^/Mpc. Observational limits placed on the 
galaxy peculiar velocities are Sv ga t axy < 600ftm.s -1 , which implies Sr = l^ft^Mpc, 
yet structures form on scales some 20 times larger than that. It appears then that 
by investigating these large scale structures we are looking directly at the primordial 
density perturbations. 

It is generally believed that the early universe was characterised by a series of 
phase transitions, during which a Higgs field ^ tended to fall towards the minima of 
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its potential. As an example, in i\<p 4 theory, the lagrangian is given by: 


£ 


1 

2 


d»<t> I 2 



<^i 2 



; 772q > 0 


(1) 


This is the familiar Mexican hat potent al. At T > T c , the critical temperature, the 
fields are in the symmetric ‘false vacuum phase with < <? >= 0. As the universe cools 
and expands through T c the <f> field rolls to the bottom of the potential developing an 
expectation value < 0 > 2 = thereby breaking the symmetry. In fact it is possible 
to show r 19j that the effective mass of the scalar field vanishes at T c with 


or 

m 2 (T ) = ml 1 - — ( 2 ) 

We will now evaluate the spatial correlation of the o field, as it determines the scale of 
fluctuations in the field 32 . The minim im work required to bring the system out of 
equilibrium for constant pressure and te nperature is the difference in the free energy, 
between the two vacuum states, AT, w th a corresponding fluctuation probability 
w cx erp(—,3AF) 

We concentrate on long wavelength fluctuations in which the d> field varies slowly 
across space. These fluctuations increase anomalously near the transition point. For 
the case of fluctuations in the symmetnc phase, < 6 >= 0,A<p = <?, then to 0(< ?*) 
we find that the change in the free energy is (from (1)) 


Af = / i rAI)|,; !+ I|^l^ (3) 

Now expanding 0(r) as a Fourier series in a volume \ 

<f> = Ek <f>ke' k - r ,4>-k = <f>l 
we find 

AF = \-J2{k 2 + m 2 (T)} ] dk \ 2 h< nee <\ 4>k j 2 > = Trrp, ( 4 ) 

2 T v W" T m ‘( T )) 

Note as T — > T c the long wavelength fluctuations increase. Writing the correlation 
function as 


G(r) =< p(r 1 )<p(r 2 ) >, r = rq - r 2 


( 5 ) 
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we use 

G(r) = 


Substituting (4) into (6) we obtain: 
T v 

G ( r ) - t -<^(- 7 ), r 

4 rrr £ 

where 




j <\ 4>k 


2 > e lfc r V 


d 3 k 

( 2 ^ 


( 6 ) 

(4 


£ 1 = m(T) 2 ^ \ A < 6 >| 


( 8 ) 


is the correlation length. ‘ This implies that domains form of size £ ~~ m -1 inside 
which 4> is correlated, but outside of which there are no correlations. However as T — ► 
T C) m — > 0. £ — * oc, it appears the domains vanish as all the fields become correlated. 
There is however an upper bound to the correlation length: from cosmology £ < 
II d 1 ( 0 ~ t , the dist ance over which microscopical forces can establish correlations 
in one Hubble expansion time. In fact there is a tighter constraint which sets the 
scale over which the domains form. As the domains form and < >— 0, there is a 

possibility that thermal fluctuations in the field could cause <f> to return to its false 
vacuum value, hence wipe out the domains. The free energy associated with such a 
fluctuation with scale £ is, using the free energy per unit volume /, 


(20 3 A/ - m~ 3 (T) 

m(T) 


X 


( 9 ) 


The fluctuation has a high probability so long as the free energy required is <C thermal 
energy available (T). The two are equal when ~ T 

1 1_ 

T 2 r c 2 ~m 2 ’ 

2 If gauge fields are present, as in superconductors, there is another relevant length which deter- 
mines the spatial correlations between the fields. It is the London penetration depth and defines the 
distance the J? field penetrates the surface of the superconductor. 8~ l — e |< (j> >j~ m^ifi). This 
scale will prove important later on. 
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(i.e. for small A, domains set in when T T c ). In fact we can then see that at thermal 
equilibrium: {( ~ [y/X < <f> >] _1 ~ [AS ’ e j _1 which is the Ginzburg length found in 
superconductivity. 

Hopefully the above arguments are convincing enough to suggest to you that it 
is worth investigating models of the ea ly universe. One such model was proposed 
by Kibble [3] in a paper that still rema ns one of the clearest in the field. At very 
early times the universe was very hot and the fields describing interactions were in 
a highly symmetric phase. However as the universe expanded and cooled, symmetry 
breaking processes would spontaneously occur, occasionally leaving behind remnants 
of the old symmetric phase, (topologica defects), possibly in the form of one dimen- 
sional strings or vortex lines [ 10 j , two dimensional domain walls or more likely three 
dimensional monopole configurations !! . In fact monopoles will always be produced 
in a GUT symmetry breaking scheme ar d this is what is commonly referred to as the 
‘monopole problem 7 , We don t see any monopoles today, however the\ would ha\e 
been produced in a great number densitv in the early universe. As they are topologi- 
cally stable, their only means of decay is hroagh annihilation or possibly gravitational 
radiation 12 . These processes are too : low to rid us of all the monopoles today, so 
how can we reconcile these apparent dis< repancies in the theories? Actually, although 
I won’t go into further details, a possible method of eliminating the monopoles with- 
out resorting to inflating them away d >es exist; inflation would also wipe out any 
cosmic strings that had been produced prior to the inflationary period, which is bad 
news for the cosmic strings. It is possible to have a sequence of symmetry breakings 
which first produce monopoles and then attaches monopole/antimonopole pairs via 
strings[13j. The monopole flux is confined to exist only on the string, and they rapidly 
come to annihilate. 

Returning to the condition for the existence of these topological defects^,, con- 
sider a gauge theory with a symmetry group G, this is the group whose elements 
leave the full potential V, invariant when acting on it. In the phase transition this 
group is broken to a sub group H as the fields pick up expectation values. H contains 
the elements of the original group G wnich when acting on the fields <^>, leave them 
with their expectation value. In fact tbe manifold of degenerate vacuum states (the 
manifold corresponding to the state of east energy in the theory), is identified with 
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the coset space: M ~ ~. 3 What then does the topology of these various coset spaces 
look like? For the answer we look to group theory. The condition for the existence of 
strings is that the first homotopy group II^ — ) = IIi(A/) be non trivial. The vacuum 
manifold M must contain non contractable loops. We will look at a specific example 
later, a more general class of models which contain non contractable loops has been 
established including some based on superstring theories[l4]. Strings are not the only 
defects that form of course. In fact if — is disconnected, then II o(M) is non trivial and 
wall like defects form where V ^ 1 ' mm inside the wall. If ^ contains non contractable 
2- spheres then n 2 (A/) 7^ 1 and the resulting defect is a monopole. 


For the rest of this lecture I will concentrate on string like discontinuities. The 
most familiar strings, i.e flux tubes in superconductors correspond to the complete 
breaking of an Abelian group G = U(1)[10-. The lagrangian for the theory is: 


c = 

4 


2 


4 ! 


/ J 2 1,1 

( <p * -7*r 


( 10 ) 


where — d M A v ^D^ — — ze.-l M ,A M is the gauge field, e the gauge coupling 

constant, A is the self coupling of the Higgs field and rj is the value of the symmetry 
breaking Higgs field. The manifold M of ground states is a circle 
M = [<£> j d> = T]t l6 , 0 < 0 < 2rr 

i.e V is minimised by <^> — tjc iS with 6 arbitrary, corresponding to the winding number 
which is an integer. We have already seen how domains could be formed during a 
phase transition. Now as the system cools below the critical temperature, (T c ~~ 77), 
the <f> field begins to fall to the minima of it’s potential. Domains form of size £ ~ 77 _1 
due to the thermal fluctuations of the (j) field. In these domains (j> points in arbitrary 
directions in M, but match smoothly at the boundary, with 6 varying so as to cause 
defects to form on the edges common to certain domains. This is easily seen. Con- 
sider one such edge where 6 varies by 2 tt in encircling the edge, i.e all around the 
edge we continuously encircle V mtn . This implies that (j) must vanish on that edge for 
it corresponds to a region where it is not in V min . Such regions line up to minimise 
the spatial gradient energy, forming a defect line or cosmic string. It corresponds to a 
thin tube of vacuum energy, V( 0 ), being stored in there. These lines where <f> — 0 are 

3 For example, if (1) G = Z 2 , (i.e V(t) invariant under <f> — * - (p ), H = 1, => V min ~ ~ ~ Z 2 

(2) G = U( 1), (i.e* — e«*), H = 1, => Vmin = £ = U{1) 

(3) G = 50(3), (Let — O fl ^ 6 ), H = 50(2), => V min = ^ = 5 a 
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either in the form of closed loops or infin tely long, for if they had ends, then it would 
be possible to move the circle (corresponding to U min ) beyond the end of the string 
and then shrink it continuously to a pe nt without having to encounter the <j> — 0 
region. Hence IIi(Af) = 1 in that region, there would then be contractable loops and 
the strings wouldn't exist 4 . 

Returning to the string solutions of (l J), we look for z- independent static solutions 
to the field equations 10b Figure (1) ‘hows the resulting solutions for i <? j and 
B_ — V A A as a function of the radial listance from the string. The width of the 
string is roughly m~ l ~ (\/A r])~ l where n 0 is the Higgs mass. The string tension, or 
mass per unit length: 



hence from U(o): 

ji ~ ( v A 7/ ) " 2 A r; 4 

- r} 2 ( 12 ) 

For example if the symmetry breaking ^cale is during the GUT era: jj ~~ 10 13 or 
10 16 Ger] then the dimensionless parameter G//, (G is Newtons constant), lies be- 
tween 10" 7 and 10 -5 . As we shall see observational constraintsT5’ place a tight 
upper bound, (not lower) on Gji of: ( 7/x < 10”°]. There is a big difference be- 

tween global and local strings. The latter, as the name implies possess a local gauge 
field whose presence results in no long : ange interactions between the Higgs fields: 
[lim r -^oo(V + ieA)<p = 0], Also the magnetic flux in such strings is quantised in units 
of [/ 5.dS = / A.dl — ~] where dS and dl are the area and line elements sur- 
rounding a portion of string. If a global symmetry is broken, there are no local gauge 
fields present, resulting in Goldstone bosons, long range forces and an infinite string 
mass per unit length. Most of the work on strings has involved local strings, apart 
from one important result we will come o soon. 

The usual method for numerically fc rming cosmic strings is referred to as the 
Kibble mechanism. A lattice of domain; is constructed typically of size — £, ( the 

4 Strings could be finite in length, by connecting them to monopoles or domain walls’ 13] 
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length scale above which the orientation of the Higgs fields are uncorrelated). In each 
domain a value of <fi £ b m ; n is randomly chosen , as this reflects (j) choosing a minimum 
energy configuration as the phase transition is passed through. After each domain 
has a value in it, look at each link on the lattice and using a prescription to smoothly 
vary the phases from one domain to the next, decide if there is a net winding number, 
anti winding number or no net winding number, hence is there a string, anti string 
or no string passing along that link. Numerical tests of this mechanism! 16] indicate 
that after the phase transition, about S09c of the string is in long ‘infinite’ string as 
long as the box in which the simulation is performed. The rest is in a scale invariant 
distribution of closed loops where the number of loops between radius r and r — dr, 

/ s dr 

n(r-)oc - 

(i.e independent of £). At high densities both the infinite and closed loops of string are 
in the form of Brownian walks of length L ~ V- More recently analytical approaches 

l '”' s 

have placed these predictions on firmer ground. First of all Mitchell and TurokT?;, 
by counting the density of states in the quantised closed bosonic string, demonstrated 
all the above results. Recently David Haws, Ray Rivers and myselfTS: using finite 
temperature field theory, have investigated the distribution of the Higgs field around 
the phase transition, also deriving the same behaviour. We are unable to derive a 
precise number for the amount of infinite string length (i.e 80%), because that is a 
phenomenon out of equilibrium where the canonical (and microcanonical) ensembles 
break down. I will spend a little while explaining our technique, although I will also 
comment on the regimes where the two results can and can not be compared. 

The simplest theory to possess vortex solutions is scalar QED 

£[4>,A] = (d^ + ieAM ] 2 i <f> | 2 | </, j 4 , m 2 > 0(13) 

The partition function is: 

Z oc J D<f>D<l>*DA(detM)ezp(-I 0 l<f>,A}) 
where det M describes the gauge fixing, and 
I 0 [<t>,A] = -f o a dTfd 3 zC E [6,A} 

with Ce the Euclidean form of the lagrangian, (i.e write <? x, r] = Y n <p n (x)exp( 2 ^ nT )) 
5 . Now to 0(A,e 2 ), for j <C 1 we find the one loop correction to the free energy, 

5 It is important to note that we are not using the mean field approach here, in that case we 
would expand about constant solutions for <f> and A 
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having integrated out the massive modes 

= -\ml I <t>o I 2 [1 “ fj\ + £ ! i 4 

where T 2 = This is then substituted back into the partition function leaving 

only the functional integration over the i ero modes remaining. 

Z = J DcpoDAo^expi -Sl A 0tl ) 

The statistical properties of strings arour d the phase transition are found by evaluat- 
ing Z with string like configurations. The dominant contributions to the integral come 

• • Q l i 

from field configurations satisfying the stationary equations S(*> 0 — 0 1 addle 

n dl \ ,, — 0 

4 8A M : saddle' M — , aiile 

d l F ;] = fi'e(c>*3,<5 - oc' ; o’) - e s .4j 6 2 


(5, - ?eA,) 2 p = -m 2 (l 


T\, \ 

■ , 


O 


o 


The contribution of any solution tc Z is found by substitution. The solution 
0 - const , .4 = 0 is a solution of minimu: l energy, so gives the maximum contribution. 
Well away from T c , Z is well represented by this term. However as T -» T c we must 
sum all maxima to the functional because the total contribution due to the large 
number of non constant field configurations becomes larger. This is easy to see, 
as J _> J c) the mass. m(T) — > 0, so it becomes possible to form strings at no 
energetic cost, and there is a second or ler phase transition. There are many string 
configurations and we must sum over them all 10, , 6 here though the solutions vary 
with temperature. 7 In particular 

lim,.-,* | 4 >(r) | — ♦ T), t ] 2 = rr, \ = m 2 ( 1 - yj) 

At the core of the string | <f> | vanishes, it’s thickness ~ m 7 1 . The magnetic field is 
determined by m" 1 , 





( 14 ) 


fl\V e should in principle sum over all the m ixima of the functional, not just the string contribu- 
tions. The approximation here has it's analogue in condensed matter with the independent vortex 
model of the A transition in liquid Helium. There the model is in quantitative agreement with 
experiment, under the same assumptions we h ive made. 

7 These solutions are strictly speaking infinite string solutions. In dealing with loops of string, if 
we are to use these solutions, then we must wo k in the regime where the radius of curvature is large 
compared to the width of the string, so the st ing is approximately straight in any given region 
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and the string has an energy per unit length a — a, -f a v 


o’. = 

o v = 


0(t] 2 (T)) scalar field 

t* 

0(rj 2 (T)) vector field 


As T Tc , cr -> 0, m 3 oc. The expression for I is substituted into Z, inter- 
actions between the strings are neglected by making them non- self intersecting, and 
placing them a distance 1 apart, where 1 is the lattice spacing which also corresponds 
to the width of a string. We then have: 

Z — £ n \V(n)exp( — 3aln) 


where \\ (n) denotes the number of configurations of a string of length nl. 

It is now possible to use standard results from polymer physics 20', to demonstrate 
the statistical properties of the strings at high and low string segment density IS . 
Restricting ourselves to non self intersecting walks, at high string densities these are 
well approximated by Brownian walks, but not at low densities. We find the partition 
function for a ‘gas' of loops is: 


Zioop — expi^Z i) 

CV x 

= exp (jp £ n ^ q ~ l e~^ nlCeff ) 

n = 1 

where C is a normalisation constant, q = | high string density, | low string density. Zj 

is the partition function of the single loops, and the effective string tension is defined 
by: 


a ~ 5, depending upon the type of lattice the cosmic strings are placed on. We 
find that in the case of large n, both the contributions from infinite string (Z x ), and 
closed loops (Zioop) diverge for: 

T > T, t 

cr l 

In( a) 

the temperature at which a e!j = 0. Above this temperature there are large fluctua- 
tions in the d> field and it is no longer appropriate to describe the field in terms of string 
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like configurations. In this region fluctuations in the energy with the canonical ensem- 
ble diverge (i.e the specific heat C v diverges as T — ► T ft ), although the microcanonical 
ensemble remains perfectly well behaved above the Hagedorn temperature[l7]. The 
reason for the breakdown in the energy is that there is an overcounting in the canoni- 
cal ensemble for the number of possible tring states becomes infinite. It is important 
to note however that this doesn’t mean the canonical ensemble always breaks down. 
For example the fluctuation in the mei n number of loops of length nl, is perfectly 
well defined all the way up to the Hagedorn temperature. In our case we find the 
fluctuations in the Higgs field are small enough just below the string temperature to 
make the canonical ensemble perfectly reasonable there. To see this, we evaluate C v , 
as this gives us the fluctuations in the energy. When this is done, we find that the 
critical exponent is j (<7 = ~) and 1 ( 7 = |), implying the fluctuations are small, 
just below T st . A very interesting resu! is obtained when we evaluate the partition 
function for an arbitrary number of d mensions d. We then have for the case of 
Brownian strings in d dimensions, q — Till 17 . Now when we evaluate C v , we find 
that the condition for it to remain finit * in the limit T — > T st is that d > 5, ( i.e for 
strings in 6 dimensions and above the fit ctuations in the energy do not diverge in the 
canonical ensemble). Of course the actu il value of T st depends also on the number of 
dimensions. 

We should think of T st as the temperature at which our strings are formed, 

T st ^ <r(T)l{T) - Wm~ l 7 - 0(1) A » e 2 

Then it follows T at < T c (i.e strings low< r the phase transition temperature) with 

1 “ = 0(A) , 771 = 777, 

= 0(e 2 ) , 777 = 777 v 

The width of the strings at formation cm now be calculated (~ mean separation) 

m,(T, t ) = 0(em,(T = 0)) , r i = m v 

= 0(v' / Am,(T = 0)) , 777 = 777, 

That is the network of strings at the pha^e transition has the separation of the centers 
of the flux tubes scaled up by a factor O ^), (e 2 < A), a result obtained by Kibblej3j, 
and discussed earlier in this lecture. T lis is a useful confirmation of the validity of 
the Kibble mechanism. What has be ome clear from the analysis is that as we 
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approach T, t , the mean field approach does not give the correct value for the critical 
temperature. We predict a second order phase transition at T, t < T c j meanjuid- As the 
symmetry is restored, overlapping strings fill the whole of space. It is possible that 
thermal fluctuations cause a first order phase transition, indeed this has recently been 
invoked[3l] in an attempt to investigate fundamental strings in the early universe. 
Above T, t ~ 0(Tg , nzburg), our model is unable to recognise string like configurations. 
As T, t is approached the system becomes dominated by infinite string, with a scale 
invariant distribution of loops approximated by Brownian trajectories. Neglecting all 
interactions, actually implies Brownian trajectories, although it would be interesting 
to discover whether the simulations of ; 33] can detect this difference in Brownian versus 
non self intersecting random walks. Below T, t the system evolves to a state with an 
exponentially suppressed distribution of large loops. 

We now move onto investigate the evolution of cosmic strings. For a string of width 
W and radius of curvature R. if R > IF then the action for a string is approximately 
a locally boosted version of the straight static solution. This is the Nambu action 21 : 

S = -fi j dadr-detg'^') ( 15 ) 

where g[ ‘J is the world sheet metric.! a.b = 0,1). In terms of the string coordinates 
A ^(<t, t) and the background spacetime metric: 


si l' = A-»A >, 


4 ) 

Vi/ 


(16) 


The Nambu action (15) also has the interpretation of being the area of the two 
dimensional worldsheet. Corrections to the action and hence to the solutions to the 
equations of motion, are of 0(^-). Since for typical cosmic strings, \\ ~ 10 -3O m, R - 
kpc the corrections are small everywhere except near a cusp where R ~ \\ [22:. 

A vital issue for string theory concerns what happens when two strings collide? 
The Nambu action breaks down here, so the full non linear field equations must 
be solved (10). This was first studied numerically for the case of global strings by 
Shellard[23], who demonstrated that strings nearly always intercommute or exchange 
partners, for relative velocities below -9. This is a very important result, as it is the 
only process by which loops can form in sufficient amounts from an initial configura- 
tion which contains nearly 80/c of its length in infinite string. Analytical approaches 
to this problem, including introducing gauge fields have been tried, but only for special 
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cases[25j. One interesting prediction that Turok and myself came up with[25], which 
proved to be the case, was in showing thf-t intercommuting always occurred when the 
strings first intersected. Whether or not they effectively reconnected later, depended 
upon the velocities of the new set of zeros that were produced in the interaction re- 
gion. Basically four zeros would be prodi ced after an initial intercommutation. These 
would then start moving out, two chasi lg behind the pair that had intercommuted 
and two in the direction the pair had b< fore intercommutation. Depending on their 
velocities the original pair would be ca lght and annihilated leaving the remaining 
pair to carry on as if no intercommutinc occurred. More often the ring of new zeros 
would not catch the original pair and would simply lose there energy as they spread 
out. The full analytical problem remains unsolved. Matzner has studied the case of 
local strings numerically and concluded they are even more likely to intercommute 
when they cross than the global case r 24 . An open question remains to be resolved- 
what happens when two superconducting strings intercommute? 


From (15), which is a purely geometrical object, we can write down the string 
equations of motion in say a flat Friedm mil Robertson Walker universe (FRW ). One 
important point to note is that the equations are independent of ji> the scale of the 
symmetry breaking. Reparameterisatio i invariance of S under a — ► a(cr, r), r —> 
t(ct,t) enables a suitable gauge to be cl osen: 

x° ~ r i\i -- 0 (17) 

where i ~ <9 r x, x = d a x. Then the equa ions of motion become, in terms of the scale 
factor a(rj) (where dt 2 = a 2 {rj)dr / 2 ), 

z + 2-(l-i 2 ) = ~(l) c = -2 -i 2 e (18) 

a e 6(7 e a 

with 


c = 



The energy in a string is 


E = 



das 


(19) 


( 20 ) 
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In general it isn’t possible to analytically solve these equations, but in certain 
regimes it can be done. 

l]. If R <C damping term, then x ~ 0, and the string is simply conformally stretched 
with a radius of curvature Hq 1 . 

2i. If R <C | = Ho\ then the string doesn’t notice the curvature of spacetime and 
acts as if it was in flat space. We can set s = 1, and the equations become 
' / 

X = X 

x.x — 0 (21) 

x 2 - x 2 = 1 

( 22 ) 


The general solution is made up of right and left moving inodes: 

x(T.a)=-a(<T — T)- L h(o r ~r} (23) 

with the gauge condition implying: 
a 2 = h 2 = 1 


in the centre of mass frame of the string. Thus c, b are closed curves on the unit 
sphere, generally intersecting if they are continuous t 26h In fact the picture for R <C t. 
is that loops of string break off the network of long string and oscillate as if in flat 
space. At an intersection x = 0, i 2 = 1 and such a point is a cusp ( defined by the 
vanishing of det^). However as x need not be continuous, such points where it isn’t 
are called kinks. Four kinks are always produced when a string intersects. 

We can summarise the early string evolution picture. Initially 80% is in infinite 
string. As the universe expands, the infinite strings chop off loops, (which sometimes 
reconnect back onto an infinite string). These loops also self intersect until a class 
of non-self intersecting loop solutions are reached. In general these loops are stable 
against perturbations which could cause them to self intersect, (i.e from say matter 
accreting around such a loop) 27]. Eventually these loops will be identified with the 
galaxies and clusters of galaxies. Many problems remain unresolved. What is the 
region of phase space which possess non-self intersecting loops? It has been shown 
to be non zero, but if we want to predict the number of daughter loops an initial 
parent loop will produce, we need to be able to answer that question. Progress 
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towards answering this issue has recently been made in that it is now possible, by 
studying just the kinks and cusps present on a loop to place a lower bound on the 
total number of times a given parent loop will self intersect [28]. How generic are 
cusps on strings, and how do we deal w th kinks are two more open questions? We 
still need to know how the gravitational backreaction from an oscillating loop affects 
its’ motion. Does it cause the loop to los; its" periodicity? As yet no exact analytical 
solution for a cosmic string in an expai ding universe has been found, so there are 
plenty of important analytical issues stil to be resolved. 

As the loops oscillate, due to their te tsion they decay primarily by emitting gra\- 
itational radiation [29]. Using the quadr lpole formula, it is possible to estimate the 
lifetime of the loop at about lU 6 oscillations. Currently work is in progress to accu- 
rately deal with the backreaction on a h-op due to it s radiation. 

The ‘loop production function’ is a:: important measurement device used in the 
evolution scenario. If the rate of loop formation from infinite strings is too small, 
then the strings quickly dominate the t< tal energy of the universe. As we ha\e seen, 
energv in strings longer than IIq 1 . ( several per horizon volume), scales as a(t), 
whereas energy in radiation scales as a ^{t). Thus ^ oc cr(t) = t. These strings 
must chop off a constant fraction of th-ir length each expansion time into loops to 
avoid this problem. The way this worts is that the correlation length, i{t) of an 
infinite string at time t, is t. This is tl e condition for the scaling solution which is 
vital if strings are not to dominate the iniverse. To see this we know at formation 

_ A 

Pitting — £ 3 £ 2 * 

Thus 

— trtn9 ~ — const if ( oc t. 

Prad 

In a scaling solution the only length sc de is Hq 1 . If events occurring 10 3:1 seconds 
after the big bang are to explain even s some 10 billion years later, we need some 
sort of scaling if we are to be able to iredict anything. Fortunately the numerical 
results do show scaling[30], and analyt cally we can explain why small loops would 
form rather than infinite string. It is due to the increased amount of phase space 
available in flat space for loops r 17,18'. 
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Another important issue which will hopefully soon be resolved in the numerical 
simulations is the role of the smallest loops that are produced. Loops decay slowly, 
so their density scales as matter, decreasing slower than in the long strings. This 
means that small loops dominate the energy density in strings. Once again the 
importance of the chopping process terminating after a finite number of intersections 
is demonstrated. We now have a distribution of loops with a scaling solution, the 
only length scale is Hq 1 and mass scale is fi. What does this model have to offer with 
regard to the formation of structure, how does matter accrete around these loops and 
how efficiently? In the next lecture I will address these issues. 


2 Lecture 2: Large Scale Structure 


Most of the mass in the universe is dark, the evidence for this comes from many 
quarters: 

1 . The fact that the spiral galaxy rotation curves stay flat much further than their 
visible radius implies ^ f> c / UJ ^ r5 - -1 - -3 

2;. Nucleosynthesis sets bounds on the baryonic matter [34], -04 < < T4. yet 

ft luminous < -02. Thus most of the baryons are dark. 

3^. Inflat ion predicts ft should be close to one. 


Fortunately particle physics provides us with many potentially exotic sources for the 


dark matter, these could be hot, (i.e neutrinos), cold (i.e axions), and thev could be 
supersymmetric 35] . 

Now the hot big bang is incomplete without a source of perturbations which were 
essential to seed structures on large scales, especially as the universe was essentially 
very isotropic and homogeneous early on. Initially we will investigate the density 
perturbations due to an infinite straight string[36]. Recall the action (15) [21]: 


s = -fi J drd(r[-detgl b ]■ 


(24) 


we can write the stress-tensor: 

2 6S 


T^(x) = -• 


g= n 


\/-g A {z)&9n*{*) 

For a straight static string along the z-axis, x(a, r) = (0,0, <r), so: 


( 25 ) 
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T M "(x) = fi8 2 (x)diag(l,0, 0,-1) (26) 

Note the role of the negative pressure term, this gives the string its’ tension and 
will prove important in the superconducting cosmic strings scenario[37j. What is the 
effect of this string on the background metric? We write, 

5V.i/(2-) = Vuv d" ^ , *U'(‘*')’ 

which is valid provided h ,< 1. In thr Harmonic gauge, (g a3 = 0), Einsteins' 
equations become 

d‘h, v = -16-G(T, u - 1 -r 1 ^) (27) 

The time independent solutions have no Newtonian potential for /too, ^33* (he V‘/i 0 o — 
V 2 /i 3 3 = 0). However we do obtain 

h n = hn 2 = SGfJiln(-) (2S) 

where r 0 is a suitable cut off introduc -d into the integral. Thus we find the line 
element becomes: 

ds 2 = di 2 - dz 2 - (1 - S Gfiln{— ) {dr 2 -■ r : do 2 (29) 

Redefining 

dr = { 1 - 4 Gfi!n(^))dr 
d(j> = (1 — 4G/i)cf^ 
to obtain 

ds 2 = dt 2 -dz 2 -{dr 2 + r 2 d p 2 ) (30) 

we see that the effect of the backreacti m of the string on the spacetime is to leave 
the metric flat everywhere except at the vertex, with an accompanying missing angle, 
0 < <j> < 2n — 87rGr/x[36] - The spacetime is conical shaped with the string passing 
through the vertex of the cone. How *ver the gravitational attraction due to the 
string on a static test particle is zero, as the particle effectively experiences a flat 
spacetime. If the string is moving relative to the test particles, they will feel the 
conical spacetime and be drawn in behind the string over an angle 2A0 = 87 rG^i, 
forming small wakes where p wa ke — 2pba k ground as matter falls in on either side of the 
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wake. I will mention observational consequences of this missing angle at the end of 
the lecture. 

\V e saw in the previous lecture how most string ends up in loops. By the time 
matter starts accreting around these loops they should be well inside the horizon, 
r <C Hq 1 , so gravity should then be weak inducing small peculiar velocities, 8v 
1. In this regime we can use linear perturbation theory [38]. Also on scales large 
compared to their size, it is a good approximation to treat the loops as Newtonian 
point particles;39j. Thus in the Newtonian approximation to Einsteins equations we 
have for a particle at r(t), (for more details see 9 ): 


r = - 

where the Newtonian potential obeys 
V 2 T = -4 ~Gp 

For cold dark matter. (CDM). mass is conserved, so in a comoving volume 
J pd\ C omoving — COTl^t 

The solution in a homogeneous background. Pe{L) only is 
inGpBZ 7 




6 

•iGpsL 


(31) 


(32) 


(33) 


Pbt — const 

In an Q = 1 universe, we find on integrating (34), 


r oc 1 a 


1 


PB 


(34) 


6 < 35 > 
which is the standard FRW matter dominated universe. Now we want to do linear 
perturbation theory about this solution. Define a(t) — (^-)f where t,- is the initial 
time. Then 


r(t) = a(t)r t (36) 

Thus x = r, is the initial comoving coordinate of the particle. Giving each particle a 
small comoving displacement v(x,t) we write 
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r(£) = a(t){x + £(x,t)} 
Mass conservation implies, 
Pid 3 Xi = p(r)d 3 r 
so using pi = pbcl 3 , this implies 

Vb 

dr 
dx 

From (37) we have: 



p(n) = PB(l-V J .v)-0(t:- 2 ) 


Defining 


it 

P 

we can write 


(p(r) - ps ) 

Pb 


-V 


v 


for the fractional density perturbation. 

To solve for r(t), hence for y(x. £). substitute (39) into (31.32) 


V 

= —i~Gp[r) 

with 

V. 

a 

and obtain 



: 4 > = 


1 L( i d '\ 

r 2 dr dr 


dQ f 

r 2 -^— = —4t tG J r 2 drpB{ 1 — flV r ,r) 


Hence we find 


V r $ = 4 ; rGp B (~ - ay) 

which when substituted in (31) gives 
H 2d ; 

p 4- — y = 4TrGp J ed’ 

~ G __ 


(37) 


(38) 


(39) 


(40) 


(41) 


(42) 
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Taking the homogeneous solution we find y « a 2 . Now the peculiar velocity of a 
test particle ( the physical velocity minus the Hubble flow) is, 

a 

v p = r r 

a 

= r — a(x + x/’) 

= ax>' 
or 

1 

V p cx 

a 

so a test particle slows down in comoving coordinates. We interpret the right hand 
side of (42) as the gravitational instability, the feedback of the perturbation on itself. 
For the case of the matter dominated FRW universe, i.e p B = the solution to 

(42) is 

v( x, t ) = A{x)t 5 — H(x)t~ l (43) 

with A and B arbitrary functions. Recalling (40) we see that — obeys the same 
equation as y*, it has the same time behaviour. Thus 

^ P / 1\ 

— oc a[t) 

P 

in the growing mode. 

What happens in the case when there is a local non-linear mass concentration? Is it 
still correct to use the above equations which are only valid in the linear regime. The 
answer turns out to be yes, provided the density perturbation is small on the surface 
surrounding the mass concentration. To see this we follow an argument presented 
in [9]. Surrounding a volume V is a comoving surface S. The average density in a 
comoving volume (for CDM) is, using (38) 

pVf = P.Vx 

= PB^Vi ( 44 ) 

where p is the average density, \\ are the final and initial volumes. Now we can 
write Vf = a 3 (V — S\ ) where 

sv = J s t- d $- ( 45 ) 
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a 3 V, 
P = PB—rr~ 
v 7 


Jcil'.dS s ^ :2s 
Pb ( 1 4 ^7 ) O V ) 


we easily see 


/ s d».dS 

7 = a u 

ex ^ in the growing mode 


(46) 


In order to apply the linear regime to n ass concentrations, we need linearity on S, 
not in V. 

Turning our attention back to peculiar velocities, recall 


v p — r - { ~)r - av. 


For the solution (43), in the growing nu.de. v p ~x t*. Now (43) also tells us that ii 
the growing mode, v oc a, whereas from 3< } we know 


, U - 

t ’ = 


Hence 


u„ = (-)> -I; = II dr 00 

p a k 

the equation used in lecture 1. The usefulness of (4i) is that in particular geometries 
Up is related to the magnitude of Del ning the Hubble velocity, u# = Hr. write 


£p _ dz 

v h r 


(48) 


In CDM scenarios we know mass is conserved in a comoving volume (33). Equation 
(46) tells us that 


8p J^-dS 

7 = vi 

For planar collapse over a length scale L. I h = HL. so 


8p _ v _ Up 

P L V H 
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For cylindrical collapse it is easy to see, (because V oc r 2 ) 

y p _ b 

v H 2 p ’ 

and for spherical collapse, ( V cx r 3 ), 

Vp _ b 

v H 3 p 

Recent observations of peculiar velocities give a unique window on ^40] . In partic- 
ular for the ‘great attractor 5 , v p — oOOkms ' _1 , at L - lOOAro 1 Mpc gives ~ — *3. 

So far we have thought only about linear collapse in a particular background. 

What happens when we include sources, (e.g loops of cosmic string)]9 ? On scales 

large compared to the size of the loop, we can treat the loops as Newtonian point 
particles. Equations (32,42) become 

= -4 ~G{pB — psou^ce) ( 49 ) 

2a : , 1 

C ^ T tnurrr (oO) 

— a ~ ~ a 

A nice trick used to solve (50) is found by recalling that we really require Using 

Gauss's theorem to rewrite (46) 

b SZ-±dv 

7 = v, (ol) 

we define S = V.v, and take the divergence of (50) 

-4-2 

& d" 2^2 ^ 2^2 ^ 4tt G£ p 30UTCe (X) t) (52) 

This has a Green's function solution, 

= J dt'G(t,t')i-6p, ourC '(x,t) 

where 

= ?(«!(*')* -r'it') 2 ). 

O 

Now for a loop with comoving trajectory Z(t,a). 

^ .. f 

Upsource\t) — p J UCT 
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which leads to the trumpet shape trajectories as the matter falls onto a loop and the 
loop shrinks in comoving coordinates 41 . 

We move on to describe briefly the ase of non - linear collapse. More complete 
treatments can be found in [9j, or in I >randenbergers review articled. Non-linear 
collapse can only be solved for special geometries. The spherical collapse model is 
ideal for the case r 3> R\ 0 o P ? where Ri oop is the radius of the loop and r is the distance 
from the loop. Then the time averaged potential from the loop is 
$ = where m is the mass of the lo«»p. For v < 1, the Newtonian approximation 

is good. In a matter dominated universe, the pressure vanishes, so one solves the 
equation for r(t ). Defining M = m — M : where is the mass inside a shell of radius 
ri, 


Mi 


2 r? 


3 ^* r ‘ 9Gt 2 


m a matter dominated universe, we find in this approximation 
r 2 GM 

— = const < U bound solution 

1 r 

> 0 escape 


For an fl = 1 universe, then, const = 0. < 11 Tie shells are bound to the loop. Imposing 
the initial condition of unperturbed Hub! !e flow . r l — If l r l — 0. the parametric solution 
for small S = ~ is 


r»( 1 — cosS ) 

r — — 

2b i 

In the early stages of collapse, 

ie = ( r Ail 1)3 

p r 


t = 


3 ti[9 — sin 6) 

At} 


(53) 


from the simple scaling of the universe as it expands. Also from M oc r 3 , 
— « ^(r ) 3 • 

p ti 


The radius r turns round when 6 = tt (r = 0), with j- = (^) 2 - 1 =4-5. 
The shell initially moves outwards with the expansion of the universe, slows down 
under gravity and eventually collapses c fter reaching a maximum radius. Roughly 
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speaking at turnaround, all the energy is potential energy, P.E 
recollapses eventually virialising (i.e kinetic energy 2^--Emoi)- 


SatlOn, T^virtal 


. Knowing this we can write 


— iutr= 8.(4 * 5). (- 7 - 


— — The shell 

T mac 

Therefore at viriali- 


(54) 


The factor of 8 comes from the collapse to and the scale factor dependence is 

from p B oc a -3 . From (54) we can estimate ~ \ vir for Abell clusters and galaxies. We 

find for the case of CDM 


— ! vir ~ 150 with 1 - Z t - 1-5 ( Abell Clusters) 

P 

— ! vtr 2= 10 4 ' 5 with 1 - Z t ~ 6-13 (Galaxies) 

P 

thus we now have an observational test for the theory(see Turok and Brandenberger 
i n 39 ). Fitting todays values for the observed overdensities in Abell clusters, ~ 
170Q- 1 , with, fl = l.h — -5, (i.e Abell clusters have only recently virialised), we 
find that for cosmic strings to be consistent, it implies Gp ~ 2.10- 6=1 . This fits 
very well with the prediction of GUT scale symmetry breaking, (p ~ 10 16 Ger), yet 
is independent of that result. With this value there are now no free parameters in 
the theory, and it is possible to show how strings predict the correct masses of the 
galaxies( see Turok and Brandenberger in 39 ). The case of Hot Dark Matter is 
different in the details of the calculation, and I advise interested readers to refer to 
[9] , and Brandenbergers article[2j for details. 

Finally in this lecture I will discuss some of the observational tests of cosmic 
string theory, tests on which the model will probably either succeed or fail. They 
all require a determination of Gp. The first test I have already mentioned is the 
gravitational lensing by strings(36], where an observer could see a double image of a 
quasar or galaxy behind a cosmic string. The images would be separated by an angle 
of 4 tt Gp ~ o'", (the missing angle), for Gp ~ 10 6 . There are five or so reported 
cases of lensing between 2.5'" - 7'". However a line of double images would be strong 
evidence for the existence of a string. 

Detection of anisotropy in the local microwave background radiation coulo be e\i- 
dence for cosmic strings, as first noted by Kaiser and Stebbins t 15.. The\ realised that 
strings produced a background temperature with steplike discontinuities on curves in 
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the sky. As with the lensing, this relies o i the canonical structure of spacetime near 
a string. Light rays travelling on either side of the string produce a Doppler shift, 
Ai/ = 87 rG/it> x in their frequency, is the transverse velocity of the string. Co- 
moving frames on either side of the string move towards each other with a velocity, 
v± ~ -4. Thus we find 


ST 

T 


= 8 ttG/u’x 

< 2 . 1 0 — 5 


which is almost at the level of accuracy ve can now observe. Recently temperature 
maps of a ‘stringy’ universe have been produced to show typical angular distributions 
of the discontinuities, hence what the t< mperature distribution may typically look 
like [42’ . 

Possibly the tightest constraint on cosmic strings comes from the gravitational 
radiation from loops of string 43.29.15 . Typically the energy density of radiation 
from a loop, in terms of the logarithmic pcctrum. 

Q g (u’) = wp(w)po\ 

(u> - R~ l is the main frequency of radiation, p 0 is the background energy density), is 

n g (u) = 2 . 10 ~ 7 h 7 o 2 . 

Now variations in the observed frequency of the millisecond pulsar, places the severest 
restriction on G/z, 

Gp < 10- 5 (^-) 8 

1 yr 

where a ~ 2tt, Tyr is the observation period in years. In principle we will soon be in 
a position to confirm or rule out cosmic trings. However in practice we are reaching 
the limits of reliability with the pulsar timing, both in our knowledge of the dynamics 
and content of the solar system and the accuracy of the atomic clocks we use to stan- 
dardise things with. Both are reaching the limits beyond which we wouldn’t trust 
conclusions inferred from their results. The latter problem could be overcome by 
using a few pulsars and comparing their relative times, as they are the most accurate 
clocks we know of. For the time being strings are still very much alive. 
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The agreement Turok found between the correlation of Abell size loops and the 
observation of the Abell clusters[44], created a great deal of interest in strings. Today 
more observations are testing the theory. Will the Abell result remain both in the 
simulations and the observations? What will the higher correlation functions look 
like. Currently a great deal of observation time is going into investigating the large 
scale structure of the universe, including things like 3pt correlation functions between 
clusters, the distribution of bubbles, voids, filaments. The next few years will be 
an exciting time in comparing the numerical simulations with observation. Work is 
currently going on to investigate the growth of matter around a network of string, 
at least in the linear regime. The full X-bodv problem is also being investigated. 
One interesting result from strings is that it predicts correlation functions should be 
universal, because of the scaling solution for strings. 

Superconducting cosmic strings also produce possible correlations between clus- 
ters of galaxies 5 . They could provide the most explosive evidence for strings, as we 
shall see in the next lecture! 

3 Lecture 3. Superconducting Cosmic Strings 

In 1984, Wit ten; 4 j realised that in some models, ordinary cosmic strings could carry 
very large currents along them, behaving like superconducting wires. The charge 
carriers accounting for the current, could either be bosonic, where a charged Higgs 
field has an expectation value in the core of the string, or fermionic where charged 
fermions are trapped as zero modes along the string. 

In this lecture I will concentrate mainly on the case of bosonic superconduct- 
ing cosmic strings, although I will discuss some of the results of the fermionic case. 
Several astrophysical consequences of the superconducting cosmic strings have been 
suggested[4,5,45,46,47j. We will review these consequences, explain how to obtain 
and understand the motion of^37,48^, and radiation from[49,48,50] superconducting 
cosmic strings. W hat will become apparent is the importance of the current to the 
equations of motion, and the fact that the total radiation from the loops is finite'50' 
contrary to previous estimates. 

First of all, let us see how superconducting cosmic strings arise[4]. In lecture 1, 
we investigated the complete breaking of an Abelian group G = U'(l) with lagrangian 
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(55) 


Instead of this £7'(1) theory, we consider : P(l)xL'(l) tm where the U'( 1 ) symmetry is 
broken in the true vacuum producing str ngs as before, but the <f> field is now coupled 
to an electrically charged x held: 


C[<t>,R».x>An] = 


-F 2 — - ft 2 

^ 1 J 1 


D,6 2 - 


D,X ; 2 -V(c5. X ) 


vhere 


U(*,x) 




(: 6 


- 2\2 


a 2 

T 


x 


-A 3 ( 6 -m*) i X 


(56) 


(57) 


with F^ u — , R^ u = d^R^ 

D^x = - ieA^x, D^4> = (< 9 „ - 13/1^ ^ 

Initially l r (l) em is unbroken, whereas f' f (l) is broken producing the usual cosmic 
strings. To ensure electromagnetism is unbroken outside the string, 


m* 


< 1 


AiA 2 m 4 

"2Af > ” 


r ZA5 ^ 

The potential is then minimised with < \ >= 0 , (electromagnetism unbroken) and 
I < (f) >| = 77 (t/'(l) spontaneously brokei ) in the vacuum. However as we know from 
figure 1 , < $ >= 0 in the core of the siring. The potential (57) for x then has the 
symmetry breaking form and forces \ 'o be non zero in the core. In fact we find 


< X > 


m in the core. It is easy to demonstrate that for some range of the 


couplings, a x condensate does exist on the string!-!, 52,51,53,54’. We will return to 
this point later. 

Our solutions are constructed arourn an arbitrary curved worldsheet with space- 
time coordinates z M (cr a ) ? where cr a ~ (r, t) are the two worldsheet coordinates. Given 
such a worldsheet we construct two spacelike normal vectors, n^(<x) (A = 1 , 2 ) which 
everywhere obeys n^x Mja = 0 . We choose them to be orthonormal (n^n^ = — Sab)- 
For any point y in spacetime closer to the string than it’s radius of curvature, we can 
associate two worldsheet coordinates cr a and two radial coordinates p A : 


y M 
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X = Xo{p) condensate. Now if Xo{p) minimises the energy of the string, so 

does e i9 Xo(p) for any real constant phase 6 , although x then carries no current. What 
Witten noticed was that at little energetic cost, the string had low energy excitations 
of the form: 


= (58) 

which led to persistent currents in the string. Substituting (58) into the static straight 
string solution along the z- axis, there is an extra contribution to the string action 
(56): 

AS X = J d 2 crd ? ; D a \ 2 

= K j dV v /3 77 a \d a B - icA a ){d b 6 + ieA b ) (59) 

where K — J d 2 p \o(p) 2 . The term ; D A Xo{p) j 2 is included in the definition of p. 
(compare (11)), and to O(p) we are left with the | D a \ \ 2 term in (59). Assuming the 
gauge field A ^ varies slowly across the string, we have used 

.4 a (<7) = <.4 M (x(cr)) 


Since 


\/- 

V ^2 


in the core of the string, and the width of the x condensate 
1 


VA 3 m 2 

from (57), we estimate K ~ j-. 
The electromagnetic current 

_ *(AS X ) 


Ja = 


6A° 


= -2 eK(d a 6 -r e.4 a ) 


enables us to rewrite (59), 

AS, = jA. I 


(60) 


(61) 
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Now the superconducting property of t ie string arises because the total change in 8 
around the loop is conserved as long as | \ '' does not go to zero, (i.e does not hop over 
the potential barrier) in the core of the string. The topologically invariant quantity 
is the winding number N, 


N = 


1 

2 x 


/ 



(62) 


where 1 is a parameter along the string Despite being the integral around a closed 
loop of the derivative of 8 , N need no ; vanish, as 8 is defined moduli 2x. N can 
be any integer. Thus from (60.62). if a situation is set up where N/0, then ^ 0 
corresponds to the ground state of the siring carrying current. The winding of 8 can 
change through the quantum tunneling >f the x field, so the current carrying state is 
metastable:4,52,51.53,54 . However for latural values of the coupling constants this 
is exponentially suppressed. 

Maxwells equations enable us to dcv rmine the gauge field .4 a in terms of d a 8. In 
the Lorentz gauge, — 0. these rear 4.50 . 

9 2 ^{y) = j d 2 ay^6' x {x[cr) - y)J"O a x u 

= J u {y) (63) 


This h as a solution in terms of the retarded Greens function, 

<?*(») = 

Z > i 

We find that for a loop with radius of curvature R, string width \Y, in the regime 
^ 1, (i.e typically), 

A»(x(cr)) ~ l-I n (I!-)J a d a x»(a) (64) 

as the leading term. The next order terms are down by 0(ln(^)) _1 . Now for a 
wire of width \V, carrying a uniformly < istributed current J, (64) tells us the gauge 
field interior to the surface .4/ = .4 0 ;'-(£)*, where A 0 is the value of the field at 

the centre of the string. So the variation of A across the string is O(^) which is 
negligible as previously assumed. Puttirg (64) into (60) we find that on the string 


J a 

A a 


K'f, 




2eh e f fd a 9 

K'ff 


7T 


)dj 


(1 


5rM#)) 


(65) 
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formula originally derived in [4j for a straight string. Typically h e jf ~ -25, ( needed 
for loops to last a cosmologically significant length of time), with a em ~ In{%) ~ 

100 . 

This still isn’t the full story to the correction to the action. So far we have obtained 
A5 X , and J 0 , .4 a in terms of d a 8. Equation (56) reminds us of the electromagnetic 
contribution of the superconducting cosmic string action. 


S tm = -j Jd'yF 

= ^ J d A yA tt d 2 A tl in the Lorentz gauge. (66) 

As previously we solve Maxwells equations for .4^ in terms of J hence J a to 
obtain’oO^ 



as the leading local contribution to the string action. Adding (67) and (61) and 
rescaling the current, we find that the action for superconducting cosmic strings in 
the absence of external magnetic fields to be 37, 


i = 


/ 


d\ 


■rl- - JaJy 


ja ~ 


J a 

Js 


( 63 ) 


Here j a is a dimensionless string current. Js will prove to be the string saturation 
current. Under what conditions is (68) valid? Certainly not in regions where R ~ W, 
i.e kinks or cusps. It ignores non- local self interactions of the string with itself, by 
radiation from one part affecting another part of the string. Actually these affects 
are down by a factor (/n(£)) _1 ~ -02. For the case of the infinite strings, we have 
dropped surface terms in S em which makes (68) invalid^55,53]. The affect is to add 
extra terms to the stress tensor. We also assumed the existence of the current on the 
string did not affect Xo(p) which was taken to be independent of J . This approximation 
is good in regions of parameter space where the Higgs coupling in (57) takes natural 
values, (i.e not small relative to e 2 ). Outside these regions, one must investigate the 
contribution case by case. Me can safely say that (68) is valid to all orders in natural 
regions of parameter space. It agrees to 0(j 2 ) with that based on a Kaluza- Klein 
construction’56 ■ , however their action is valid only for small j. 


30 



ORIGINAL F AGE IS 
OF POOR QUALITY 


Before discussing solutions to (68), ve will review some of the proposed astro- 
physics! consequences of superconducting cosmic strings. Chudnovsky et al.[45] con- 
sidered their affect on plasma, estimati lg the synchrotron emission and radiation 
from shock heating. Massive strings w.th Gji ~ 10 6 , have been invoked by Os- 
triker, Thompson and Witten, [OTWj'o', as the power sources for giant explosions. 
We will look more closely at their very ittractive scenario here. The recent survey 
of the northern sky[8] seems to indicate that galaxies are found on the surfaces of 
bubbles of the size~ 20h^Mpc, the surfaces surrounding voids which are apparently 
empty of galaxies. So far there has beei no real explanations from ordinary cosmic 
strings for the existence of these bubble, or voids, (I should point out here, as more 
data is being taken, some of the voids ; re beginning to close in, proving to be not 
as empty as was first thought). Howevei superconducting cosmic strings do seem to 
offer a natural explanation for these voh s. We know that as loops oscillate, g,a\ita- 
tional waves are radiated with a luminosity, L g = - G /*' * • causing the loop to shrink. 


(-j g ~ 50 - 100 independent of the loop size). Now because superconducting cosmic 
strings carry an electrical current, J < - e v R- an oscillating loop will also 

emit electromagnetic radiation. We can renerally write the electromagnetic luminos- 
ity, L tm = 7 em jVi where j = ' sm - 10. again independent of size. Later 

we will see how the string equations of motion modify the power radiated. OTW 
show that for Gp ~ 10‘ 6 , L em > L g if j - -2. We have already seen from (62,65), 
once a loop is initially threaded by an ixternal electromagnetic field which is then 
removed, a current is induced in the loop, trapping a fraction of the flux of order 
one. Conservation of N ensures that flu : is conserved. Now as the loop oscillates, it 


shrinks; as the radius Ri of the loop dei reases (62) tells us J will scale as J oc R l . 
Thus in theory even if the initial current <C */ maT ,Z em eventually dominates L g . 
The very low frequency waves emitted by the loop can not propagate through the 
plasma, the ambient plasma frequency s 13 or so orders of magnitude greater than 
that of the emitted waves. Thus each loop blows a bubble around it to a final radius 
R. « ii M^(l + z d )^Mpc{ 5 ]. The loops heat up their surroundings, generating large 
dense shells of gas, the galaxies form cn these gravitationally unstable sheLs with 
their present distribution representing ‘he distribution of bubbles when fragmenta- 
tion occurred. The technical details of t le calculation can be found in ;5>. They find 


that for 
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R, ~ 10 — 20 h~ 1 Mpc, in agreement with [8]. Dark matter remains inside the voids, 
if it is cold then it can not catch up the expanding shells of gas, indicating that 
local dynamical measurements of Q in the vicinity of galaxies and clusters will be a 
significant underestimate of the true value. As the loop continues to shrink, J —* 
J max • The string saturates, producing high energy cosmic rays, charged particles, 
photons, neutrinos and possibly as yet unknown physics. Bose carriers lead to a first 
order phase transition and a non superconducting state, at J mQI ~ 10 2 ° amperes 
for GUT models. In a typical case the energy released electromagnetic radiation 
is 10 49 ergs * 1 or 10 66 erg in total. When the current saturates, a loop will emit 
vast amounts of particles. The loop may be seen as a high energy x-ray source at 
z 10 — 50 . Such loops will also contribute to the hard x-ray and 7-ray backgrounds, 
aswel as the 10 2p ev cosmic rays. One possible problem this scenario faces, is how to 
obtain the primordial magnetic fields. With Gg ~ 10 6 , they require the primordial 
energy density in magnetic fields ~~ 3 . 10 ~ 9 of the radiation energy density to induce 
large enough currents. No natural mechanism exists to generate such large magnetic 
fields. Random fluctuations in the x field on the string, will induce a very small initial 
current as we shall see. 


\ ilenkin and Field [ 46 j proposed that loops of superconducting cosmic strings which 
possess cusps emit short bursts of highly directed electromagnetic radiation from these 
regions, possibly accounting for the jets observed in quasars, formed from accelerated 
particles. The loops of cosmic strings provide the central engines of the quasars. The 
calculation is based on the assumption that the mean power radiated from super- 
conducting cosmic strings at current j is, P oc j$, and is dominated by short cusp 
bursts[ 48 , 49 ]. This relationship and the effect of current on cusps will be discussed 
later. The total pow r er emitted is estimated at ~ 10 62 erg. 

Returning to the action for the superconducting cosmic strings, (68), we will now 
look at the equation of motion and obtain solutions[ 37 ]. Then we will look at the radi- 
ation from these solutions and compare the answers with previous estimates! 48 , 49 , 501 . 
The variation with respect to 6 and yield 


d a {^l ab j b ) = 0 S a (v^( 7 ai + Q ah )d b x ") = 0 


where the world sheet energy momentum tensor is defined by 


( 69 ) 
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(70) 


Choosing the orthonormal gauge, 7<r T = 5, 7<rcr + 7 tt — 0 then using (65) in (69) we 
obtain 


(d*. 02)0 = 0 or 0 = f{<r-T)-g{a-T) (71) 

with f and g arbitrary functions. Xow tie orthonormal gauge is invariant under the 
coordinate transformation, 

2 

a - t -» a -r f = —f{( 7 - r) 
cr — r—+a — T = -- g{<r - r) 

so in those coordinates, 6 = Aa only, p ovided the current is non-zero everywhere. 
We have from (65) 


jo = 0, 



(72) 


and (69) becomes 

d T [( 1 -7 n j 3 )drX M ) = da{ 1 -7 u j 2 )0 o x- 


(73) 


7 


11 




2 ) -1 


(74) 


Now if (c^z) 2 + j 2 — 0, then we have the solution x° = jr, z = z(<x), an arbitrary 
function subject to the above condition . These arbitrary static curves are called 
‘springs’[37]. What is happening is th? t the positive pressure contributed by the 
current cancels the string tension. This nay be seen from the stress tensor 


Srjuu 

= n J dV v '7IX( 7 Qi> + 0 ob )0 .x (j a i ,x^ 4 (x^ - x 4 (cr )) (75) 

In our gauge, for a loop of large radius < f curvature, we approximate it as a straight 
line in the z- direction, i.e x 4 = (r, 0,0, r). We then have 

T 4 " = gS 2 {x)diag{\ - j 2 , 0,0,1 - : 2 ) (76) 
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This result differs from[55], as mentioned earlier. Really there will be surface terms 
contributing to T rr , Tee in the case of the infinite string. For j = = 1, the 

negative pressure or string tension term is cancelled, leaving effectively a line of 
pressureless dust. This is the significance of J5. It appears that as the current 
approaches the critical value, springs will form, they behave as matter and will not 
decay for they have stopped oscillating. This is potentially disasterous and in fact 
was first mentioned by OTW[5,37b If the current density is non zero, we can not 
have ‘cusps’ in superconducting cosmic strings, for as x — > 0,7 n — > oc, yet 7 11 j 2 = 1 
is a gauge condition. 

Another class of solutions are oscillating loops which are also solutions to (73)[50 . 
For x 2 = const , Xq = r, we have. 


where. 


v 2 = 


1 - ^ 


with the general solution 


1 . 


= UM 0- ~ l ' r ) - ~ vr Y: 


The gauge conditions imply 


a 2 = i = C 2 = 2 


1 + V 


a 


( 77 ) 

(78) 


when we impose the condition d.b = 0. Thus a and b are curves lying on a sphere of 
radius (. They are periodic with zero centre of mass. Solutions of this type have been 
found for ordinary cosmic strings [49]. They possess kinks, where the x derivative is 
discontinuous. The simplest kinky loop has a being two points, the north and south 
pole, b two antinodal points on the equator. From (77,78) we find x 2 = , hence 

we have 

2 ^(l + 8P)-(l + 2p) 

v — 

2j 2 
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The effect of the current has been to reduce the velocity of wave propagation on the 
string to v 2 < 1. As j increases, v decreases until v = 0 when j 2 = 1, the kinks do not 
propagate around the loop and we are ii the realm of springs. The fact the current 
slows down the kinks is enough to make he power radiated from the oscillating loop 
finite. Previous estimates at the electromagnetic radiation from loops have assumed 
they obey a Nambu trajectory. This has led to infinite answers for the total power, 
hence the introduction of cut offs into the answers 49,48b Yet we have seen the 
dramatic effect the current has on the m >tion of the loop. Using the solution (77,78) 
we were able to evaluate exactly the power radiated from a kinky loop, and obtained 


dP 

In 


ocj 2 (l -j 2 yin(-:) 


which is finite for both small and large j. In fact the total power radiated is also 
finite, tending to zero as j — > 0, and as j — > 1. The explanation is simple. It comes 
from the fact that kinks slow down pro iucir.g springs as j — > 1. The overall result 
turns out to be far less power than pre ious estimates gave, making the quasar jet 
scenario difficult to reconcile with these string solutions. 

We have seen a simple argument that radiation from loops and their consequent 
shrinkage leads to a build up in current. Basically because the total winding number 
of the 6 field is fixed, (barring tunneling events), j oc Y$. As tne loop loses energy 
it shrinks, therefore Y 8 rises. This resu t is important in the OTW scenario, where 
whatever the primordial magnetic field [ B primordial) is, eventually a loop will radiate 
away lots of electromagnetic radiation. This power, P em ~ FcmfJ-j 2 dominates the 
gravitational radiation, P g ~ T s G/x 2 for j > x'Gfi ~ 10 -3 . Thus j > 10 3 for any 
interesting effects. However when j = l, springs form, so we require 10 -3 < j e < 1. We 
know that as j — * 1 the Nambu action acquires corrections, what about the region 
close to the lower limit? Recall J a = 2 eA e ffd a d. If any region has V6 > Y0 C , it will 
turn critical and lose its’ current. This translates into if | x |< in that region, 
current will be lost. Now for small j c the motion of the loop should be nearly Nambu, 
allowing cusp like regions where | x ~ 0 These regions will lose current. Meanwhile 
the loop shrinks an amount A L ~~ TgGfiL every period, L is the length of the loop. 
Thus there is a current gain due to radiation from the loop and a competition 


d[ t __P_ 

it j,L 
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dj j j 3 

Jt lflairi= TGfl L +Tem J 


There exists a stable fixed point at j ~ lO -4 ^ when ^ = 0 [50] . This value is too 
small to cause the explosions required, hence would imply large initial currents or 
primordial magnetic fields, as loop shrinkage appears to lead to current loss, not 
buildup. 

Many of the possible problems with superconducting cosmic strings lie in the 
observation that they like to form springs. This of course is model dependent, so 
given the region of parameter space in which superconducting cosmic strings form, 
what fraction of these solutions will be springs? Many authors have tackled this 
issue L 52,53,51 ,54] , and come up with a wide cross section of conclusions. A common 
result is that springs can form, but do not have to. However, it appears that if 
we demand that superconducting cosmic strings remain around long enough to be 
cosmologically interesting (A^ > 20), this constrains the parameters in the Higgs 
potential, the result being that such loops with long lived currents are generically 
springs (i.e exist for a wide range of 

Recently solutions have been discussed in which the string possesses charge jo 
aswel as current jij57j. These objects are a kind of cross breed between the usual 
spring solutions and the noil-topological soliton solutions, in that the presence of 
a conserved charge enables the strings to be stabilised for a much lower current 
than is required to form springs (for couplings A 3 1). These Xorton’ solutions if 
stable against decay, would behave as matter, and very quickly come to dominate 
the energy density of the universe, a big problem for the explosion scenario, and for 
superconducting cosmic strings in general. Again by carefully choosing the coupling 
constants, like springs, vortons could be avoided. However cries of ‘fine tuning’ may 
well then be heard. 
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